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If the phase transition of QCD at nonzero temperature is 
dominated by the (approximate) restoration of chiral symme- 
try, then the transition might be characterized using a gauged 
linear sigma model. Assuming that vector meson dominance 
holds, such sigma models predict that at the temperature of 
chiral restoration, the pole mass of the thermal p meson is 
greater than that at zero temperature; in the chiral limit and 
in weak coupling this mass is ~ 962 MeV . The width of the 
thermal p — a\ peak is estimated to be about 200 — 250 MeV . 
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When the quark-gluon plasma is produced by the col- 
lision of large nuclei at ultrarelativistic energies, such as 
at RHIC and LHC, the crucial question is how to detect 
its presence as the plasma expands and cools into ordi- 
nary hadronic matter. A promising signal is to look at 
the production of dilcptons, since they escape from the 
fireball essentially without interaction. The most promi- 
nent feature of the dilepton spectra are the peaks from 
their coupling to vector mesons. 

Vector mesons can be classified into two types. For 
mesons such as the p, their lifetime is so short that they 
decay within the plasma. Consequently, the shift in their 
mass and width from interactions in the plasma — the 
nature of the "thermal p" — are in principle observable 
H . The second type of mesons are those whose lifetime 
is so long that they decay outside of the plasma, like 
the J/^>. Then any shift in the mass or width is not 
observable, but one can measure a relative depletion in 
the height of the peak ]!| . 

In this Letter I investigate the nature of the thermal 
p within the context of a gauged linear sigma model || . 
Several other authors have conducted similar studies in 
these [|J, H and other ||- jl^] models. The principal 
point herein is that at least in weak coupling, a general 
feature of gauged linear sigma models is that at the point 
where chiral symmetry is restored, the mass of the ther- 
mal p is greater than that at zero temperature. The shift 
in the p mass can be relatively large, on the order of T x , 
where T x is the temperature for the restoration of chiral 
symmetry. 

I work with two flavors, assuming that the effects of 
the axial anomaly are always large, so the global chiral 
symmetry is SU(2)i x SU(2) r . Introducing the matrices 
t° = 1/2 and t a , tr(t a t b ) = S ab /2, the scalar field $ is 

$ = at + iir-t ; 



7? is the J p — isotriplet pion field and a a + isosinglet 
field. For the left and right handed vector fields I take 



where lu and p are 1~ fields, and /i and eti are 1 + fields. 
According to the principle of vector meson dominance 
||, the dimensionless couplings of the vector fields to 
themselves and to $ are exclusively those which follow 
by promoting the global chiral symmetry to a local sym- 
metry. Introducing the coupling constant g for vector 
meson dominance, the appropriate covariant derivative 
and field strengths are L>^$ = - ig{Af $ - <J>A£?) 
and Fff = d»A\ r - d v A? r - ig[A? >r , A\ r \. The effective 
lagrangian is then 



C = tr (JZ^$| 2 - ^i 2 |$| 2 + A(|$| 2 ) 2 - 2fct°$ 

+\(Fr) 2 + \{Fn 2 +™ 2 m? + {A^ . a) 

Including g, the parameters of the model are a mass 
squared — p 2 , which drives spontaneous symmetry break- 
ing at zero temperature, a dimensionless scalar coupling 
A, a background field h to make the pions massive, and 
a mass term ~ m 2 for the gauge fields |b|] . Much of the 
physics of this lagrangian can be understood from the 
kinetic term for the scalar field. 
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f(l^i 2 ) = 2r (7+ "" i 



+ (c^tt + gp^xn- gafaf + g 2 (a 2 + vr 2 ) (/f ) 2 ) (2) 

Because it couples to the (isosinglet) current for fermion 
number, the w M field drops completely out of (||). There 
are interactions of ui^ due to effects of the anomaly, but 
these are neglected in this work. 

I stress how remarkable the principle of vector me- 
son dominance is. If one constructs the most general 
lagragian consonant with the global chiral symmetry of 
SU(2)i x SU(2) r , then instead of a one coupling con- 
stant g, many more dimensionless coupling constants are 
required ]l4| . Vector meson dominance limits the break- 
ing of the local chiral symmetry solely to soft mass terms 
p3[ , such as that ~ m 2 in (Q) . As I discuss at the end of 
this Letter, if the principle of vector meson dominance is 
abandoned, then very different predictions follow. 
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Of course the price paid is that the theory is not per- 
turbatively renormalizable. For a vector field with mass 
to, in momentum space the propagator is A flu (P) = 



(S^ -P»P v /m 2 )/(P 2 



TO 



which is <~ 1 and so badly 



behaved at large P. In the present analysis this lack of 
renormalizability is inconsequential. This is because I 
assume that I am always in a regime where the temper- 
ature T < T x <c to, and for such low temperatures the 
effects of quantum vector fields should be temperature 
independent. 

When spontaneous symmetry breaking occurs, so a — > 
(7q + Cj the vector meson masses are 



(3) 



Further, from (|2j) Co 7^ generates a mixing between 
the field with c? M 7r. This produces a type of "par- 
tial" Higgs effect, whereby the standard results in a linear 
sigma model are modified by ratios of m ai /m p : 



-co, m, 



m ai h 2 h 2 

rn a = — + 2Xa . (4) 



to 2 , cr 



00 



In MeV I use the values /„■ = 93, to^ = 137, to p = 
770, and m ai — 1260. Notice that the value of ratio 
Tn ai /m p ~ 1.6 is significantly larger than one. These 
values determine 00 = 152 MeV, g — 6.55, h = 
(102 MeV) 3 , and m = 770 MeV. The values of the re- 
maining parameters depend upon the value of m a . I 
choose two representative values Q; = 600 MeV 
gives A = 7.62 and /z = 412 A/eV, while to^ = 1000 MeV 
gives A = 21.4 and /ti = 700 MeV. With these values 
of A and g the theory is manifestly in a strong coupling 
regime. Nevertheless, to gain a qualitative understanding 
of the physics I work to lowest order in a loop expansion. 

In weak coupling it is easy to compute the thermal 
masses at the temperature of chiral symmetry restora- 
tion, T x . For simplicity I work in the chiral limit, h — 0, 
where T 2 = 2a 2 , so T x = 215 MeV |§, @- At T x 
I can compute in the symmetric phase, working from 
above. A technical but crucial point is that it is neces- 
sary to compute the self energies not at zero momentum, 
but on the relevant mass shell, since this is what deter- 
mines the coupling to dileptons. Consequently, instead 
of the low momentum limit of the self energies, one is 
interested in their limit for large momentum P 3> T. 
Calculation shows that the p and a\ self energies are 
each IP" = (5^ - P^P V / ' P 2 ){g 2 T 2 /6), while the fx self 
energy is IP 1 ' = 5^{g 2 T 2 /3) at large P > T. Using 
T 2 = 2a 2 , and (|J), in weak coupling at the critical tem- 
perature the pole masses in the vector meson propagators 
are given by 

m 2 p (T x )=m 2 ai (T x ) = 1 (2m 2 p + m 2 ai ) = (962 MeV) 2 , 
m 2 h (T x ) = \{rn 2 p + 2r<) = (1120 MeV) 2 . (5) 



On the right hand side of (||) and henceforth, whenever I 
write a mass such as m p or m ai , implicitly I am referring 
to their values at zero temperature; any thermal pole 
mass is denoted by m p (T), etc. 

Since in (|J) the lo field docs not interact with the scalar 
fields, the w mass does not move, to 2 (T) = mj; rrv^iT) 
only shifts from effects of the anomaly. At the very least, 
it is apparent that the near degeneracy between the zero 
temperature masses of the u> and the p, and the ai and 
the /1, is badly broken at nonzero temperature. 

The width of the p can be computed by standard means 
p8[ ; at one loop order the only available mode is p — > -kit. 
For a p decaying at rest, 



r x = 
p 



48tt 



2n(TO*/2)) u p> 



4 ( TO X)2)3/2 



(mp) 2 



(6) 



Here m£ = to^(T x ) and m£ = m p {T x ) are the thermal 
pole masses at T = T x , and = T p (T x ). This is just the 
standard formula for the decay width of the p, except that 
there is a factor involving the Bose-Einstein distribution 
function, n(E) = l/(exp(E/T) — l), from stimulated pion 
emission in a thermal bath. At zero temperature, (|^) 
gives a decay width that is about 20% too large, r p (0) ~ 
179 MeV instead of the experimental value of 150 MeV. 

To obtain a somewhat realistic estimate of the width 
of the thermal p, the nonzero mass of the pion must be 
included. The full problem with h ^ and T ^ is 
rather complicated, since toJ ~ T x . I adopt an approx- 
imate solution: the thermal effects are computed in the 
high temperature limit, including only the terms 5C = 
(XT 2 /2)tr{\^\ 2 ) + (g 2 T 2 /12)((p^ 2 + (ai) 2 ). When h^0 
the definition of T x is ambiguous; I define T x as the point 
where m a {T) has a minimum with respect to T. Doing 
so, for m a = 600 MeV I find T x = 226 MeV; at T = T x , 
fx = 32 MeV, = 978 MeV, = 1002 MeV, 

ml = 185 MeV, ml = 221 MeV, and VI = 278 MeV. 



221 MeV, and 
1000 MeV I find: T v = 221 MeV; at T = T, 



For m a 

fx = 23 MeV, mx 
toJ 



971 MeV, m* 



983 MeV, 

217 MeV, ml = 263 MeV, and = 248 MeV. If 
I assume that the p width is too high by the same amount 
at T x as at T = 0, and so should be corrected by a factor 
of 150/179, I obtain = 233 MeV for to ct = 600 MeV 
and r£ = 208 MeV for m a = 1000 MeV. 

The form in which I have written (|^) is a bit mis- 
leading, in that at leading order in weak coupling I can 
eliminate g entirely, to write expressions for the masses 
at T x solely in terms of the zero temperature masses. 
Nevertheless, it should be emphasized that this is a trick 
only of results to lowest order; the corrections to @ and 
(H) are a power series in g 2 and A, and so large. Thus 
the above numerical values are not meant to be taken as 
predictions, but only as suggestions of the magnitude of 
the possible effect. Perhaps, however, the qualitative fea- 
tures of a weak coupling analysis are reasonable. At zero 
temperature the splitting between the p and ai masses 
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are driven entirely by spontaneous symmetry breaking; 
it is sensible that the thermal fluctuations which restore 
the symmetry are of the same order as the shift upward 
in the (thermal) p mass. Similarly, while thermal broad- 
ening can be very significant if the p mass decreases, if 
the p mass increases these effects are naturally small, 
since then the 7r's are energetic, with momenta signifi- 
cantly larger than the temperature. One effect which 1 
have neglected which increases is the thermal width of 
the 7r's; however, a more realistic value of T x is probably 



lower than the above |16|, which lowers 



It is also of interest to compute the shift in the pole 
masses at low temperature. In the chiral limit we can 
make comparison with a general analysis of Eletsky and 
Ioffe , who show that the shift in the pole masses van- 
ishes to order ~ T 2 about T = 0. In gauged sigma models 
this holds for both the p and a\ masses || . The first non- 
leading terms in the pole masses for the transverse fields 
are, in the chiral limit, 



a 9 Vr 4 K,( 3m +V) " 

*pv- i - m P ~ ' * I — — —~ 3 



m 2 JT) 



45m 2 



{m 2 ai - m 2 ) 2 



m 2 ai (T) 



g 2 n 2 T 4 /4<(3<+V) 



45m 2 



( m2 ai 



m 2 ) 2 



2m 4 




(7) 



where p 2 is the spatial momentum squared of the field. 
That is, while by the time of the chiral transition the 
thermal p mass goes up, and the a\ mass down, about 
zero temperature they start out in the opposite direction: 
the p mass goes down, and the a\ up! 

Putting in the values of m p , m ai and g, at 
zero momentum, p = 0, I find that (m 2 (T) — 
m 2 p )/m 2 = ~{2MT/m p )\ while (r<(T) -m 2 ai )/m 2 ai = 
+(3.16 T/m p ) 4 when m a = 600 MeV, and {m 2 ai (T) - 
m 2 ai )/m 2 ai = +(3.17 T/m p ) 4 for m a = 1000 MeV. These 
values are interesting because the coefficients of T/m p on 
the right hand side are relatively large: if we push them 
well beyond their range of validity, to T ~ 200 MeV, they 
suggest that the shifts in the thermal p and a\ masses can 
be significant, on the order of T x , as found in (||). 

The shift in the thermal masses at low temperature 
can also be computed away from the chiral limit. When 
TOtt / I find that the p mass does not shift to ~ T 2 , 
but the a% mass does, 



4co 



9 2 mlT 2 
Am 2 



+ 



(8) 



As for the ~ T 4 term in the chiral limit, (Q), when 
m v + 1 the a\ mass starts out by going up at low 
temperature. In QCD, except at the very lowest tem- 
peratures, this correction is small relative to that in 



(@): - m 2 J/m 2 ai = +(.46 T/m p ) 2 for m a = 

600Af e y, and {m 2 ai (T) - m 2 ai ) /m 2 ai = +(.27 T/m p ) 2 for 
m CT = 1000Me7. 

I conclude by discussing the relationship with other 
approaches. By using a gauged linear sigma model for 
T < T x , implicitly I am assuming that the behavior 
of QCD at nonzero temperature is dominated by the 
restoration of chiral symmetry, and not by deconfinc- 
ment. This accords with current numerical simulations 
of lattice gauge theory |l7]] , which indicates while there is 
no true phase transition in QCD, it lies close to a chiral 
critical point |lf| . 

In contrast, if the phase transition were dominated by 
deconfinement, then as argued initially in fil, it is con- 
ceivable that the thermal p mass decreases with increas- 
ing temperature. For example, sum rule analyses of the 
phase transition can be construed as dominated by de- 
confinement. Generally, such analyses find that the ther- 
mal p mass goes down as T goes up 0] (see, however, 
[jlO|| ) ; about zero temperature, ref. JO] find that both the 
p and a\ masses decrease to ~ T 4 , contrary to (0). Using 
the experimental phase shifts, Shuryak and Thorsson jf| 
also find that the p mass decreases, by a small amount, 
at T - T x . 

Following Georgi, Brown and Rho, and others ||, have 
analyzed a sigma model where the p mass decreases 
monotonically with temperature. While their analysis 
uses a nonlinear sigma model, it can be reexpressed in 
terms of a linear sigma model. Assume that the explicit 
mass term for the gauge fields ~ m 2 in (|l]) vanishes, and 
that instead the local chiral symmetry is broken only by 
the term such as C K = K£r(|$| 2 )tr((^) 2 + (/l{?) 2 ), where 
k is a dimensionless coupling constant. With such a term, 
up to T x the p mass does decrease uniformly with temper- 
ature; an easy calculation shows that in the chiral limit, 
m 2 p {T x ) = m 2 ai (T x ) = (2/3)m 2 = (629 MeV) 2 . However, 
setting m = and including C K manifestly violates the 
assumption of vector meson dominance, since then the 
local chiral symmetry is broken by a term with a dimen- 
sionless, instead of a dimensional, coupling constant. 

In other words, which way the thermal p goes depends 
crucially upon whether or not vector meson dominance 
applies at nonzero temperature. If vector meson domi- 
nance holds, the thermal p mass goes up by T x ; without 
vector meson dominance, there is no unique prediction. 

To understand the relationship to the theory at T > 
T x , it is necessary to remember that a constant feature 
of the lattice results O] is that independent of the or- 
der of the phase transition, uniformly there appears to 
be a large increase in the entropy in a narrow region of 
temperature. Such a large increase in entropy cannot 
be described by the kind of gauged linear sigma models 
which I have been using. Consequently, I presume that 
such sigma models are valid only to a temperature just 
below T x , but not above. 

In heavy ion collisions at ultrarelativistic energies, 
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then, if a mixed phase lives for a long time and dominates 
total dilepton production, a two state signal should ap- 
pear in dilepton production. From the quark-gluon phase 
at T — T+, dilepton production is dominated by the 
quark quasiparticles ^0|, presumably concentrated in a 
region below the zero temperature p peak. The hadronic 
phase at T — T~ generates a thermal p peak; the posi- 
tion of this peak is model dependent, lying either above 
or below the zero temperature p peak, depending upon 
whether the assumptions of Q, Q, and this work, or 
those of §, §, |],and @, apply. 

Whichever scenario applies, theoretically there are nu- 
merous indications that if it is possible to resolve rela- 
tively wide structure in dilepton production in ultrarel- 
ativistic heavy ion collisions — on the order of T x ~ 
200 MeV — then it might well reveal novel structure. 
While experimentally this is an extremely difficult task, 
the possible rewards appear well worth the effort. 

I happily (if belatedly) acknowledge that an inspira- 
tional colloquium on the quark-gluon plasma by W. J. 
Willis at Yale University in 1981 originally jjj stimulated 
my interest in this problem. During the present inves- 
tigation I benefited from discussions with J. Bijnens, V. 
Eletsky, T. Hatsuda, S.-H. Lee, M. Rho, E. Shuryak, A. 
Sirlin, C. Song, L. Trueman, A. Weldon, and especially 
S. Gavin. This work is supported by a DOE grant at 
Brookhaven National Laboratory, DE-AC02-76CH00016. 
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